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ABSTRACT 

Measurements of open charm hadro-production from CERN and Fermilab 
experiments are reviewed, with particular emphasis on the absolute cross 
sections and on their A and y/s dependences. Differential pT and xp cross 
sections calculated with the Pythia event generator are found to be in rea- 
sonable agreement with recent data. The calculations are scaled to nucleus- 
nucleus collisions and the expected lepton pair yield is deduced. The charm 
contribution to the low mass dilepton continuum observed by the CERES 
experiment is found to be negligible. In particular, it is shown that the 
observed low mass dilepton excess in S-Au collisions cannot be explained by 
charm enhancement. 
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1 Introduction 



In order to understand whether the data cohected in nucleus-nucleus col- 
lisions signal new physics such as quark-gluon plasma formation, there are 
essentially two possible procedures. The first one consists in using Monte 
Carlo event generators to calculate what should happen in conventional sce- 
narios. A significant difference between prediction and experiment would 
indicate new physics. The standard event generator codes, however, simu- 
late nucleus-nucleus collisions only in what concerns soft processes. Physics 
results based on conventional hard processes in nucleus-nucleus collisions, 
on the other hand, can be predicted by extrapolation from pp and pA data. 
This second method has been used in the study of high mass dimuon pro- 
duction. The present understanding of J/V', ip' and Drell-Yan production in 
pA interactions, from data collected at CERN and at Fermilab, allows to 
estimate the behavior expected in S and Pb induced collisions 

Open charm production is a hard process and the presently available 
nucleus-nucleus event generators cannot be used to estimate production 
yields or contributions to the dilepton spectra measured at SPS energies. 
The SPS heavy ion collaborations NA38 and HELIOS-3 have used the 
Pythia Monte Carlo event generator Q to calculate charm production on 
the pp level and used a scaling which is linear with the nuclear mass num- 
ber A to compare the estimated values to the data collected with sulphur 
beams [^,^. An initial attempt to assess the charm decay contribution to 
dilepton spectra at higher masses has been presented by A. Shor ||5|. In this 
paper we use Pythia to evaluate whether the low mass electron pair excess 
observed by CERES might be explained by charm enhancement in heavy 
ion collisions. 

The theoretical understanding of the features seen in the charm data for 
pp and pA collisions has made substantial progress in recent years but is 
still facing some difficulties (see Q for a recent review). The availability 
of NLO calculations has not eliminated the need for uncomfortably large 
values for the intrinsic transverse momentum of the incoming partons 
(the '/ct kick'). Adding together the uncertainty in c-quark mass {rric rang- 
ing from 1.2 to 1.8 GeV/c^), values of {k^) that can go up to 2 GeV^/c^, 
doubts on which fragmentation scheme better describes the hadronization 
process, and the ignorance of higher order effects, especially at low trans- 
verse momentum, we get the feeling that the theoretical understanding is 
still somewhat unclear. For a detailed discussion of these theoretical uncer- 
tainties including a complete list of references see Q . It is not our intention 
to clarify this situation. Rather we wish to evaluate whether a Pythia 
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based model is able to describe reasonably well the available data on open 
charm production. If so, we can use it as a good basis to predict, in the 
absence of new physics, charm production in nucleus-nucleus collisions and 
to estimate the charm contribution to dilepton spectra. 

In Section 2 we compare the results provided by Pythia with the recent 
FNAL and CERN data on open charm hadro-production, with particular 
emphasis on the energy dependence of the charm cross section, on the dis- 
tributions in transverse momentum px and the Feynman x variable xp and 
on the correlations between charmed hadrons. We then compute the lepton 
pair yield from simultaneous semileptonic decays of D and D mesons and 
Ac baryons and compare it to the experimentally detected inclusive dilep- 
ton yield in p-p and p-A collisions. Scaling the results to nucleus- nucleus 
collisions and comparing them to data as well as the possibility of enhanced 
charm production in these collisions are discussed in Section 4. 

2 Charm production in pp and pA collisions 

The recent Fermilab and CERN experiments on the production of D mesons 
in proton and pion induced collisions provide a rather detailed picture of 
charm hadro-production. Fig. |I| collects the forward (xp > 0) D meson 
production cross sections as measured by several experiments, as a function 
of ^/s, the center of mass energy. Where appropriate, the values have been 
adjusted to reflect the current knowledge of branching ratios Q. The cross 
section values for proton induced reactions, published by experiments NA27, 
E743 and E653 for all xp, were divided by 2. The results obtained with nu- 
clear targets were divided by the target atomic mass number A. The error 
bars reflect statistical and systematical uncertainties added in quadrature 
and include the uncertainty on the branching ratios. The order of magnitude 
discrepancies found in earlier measurements have disappeared and the en- 
ergy dependence of the data is rather consistent with PYTHIA based QCD 
predictions (solid lines in Fig. |l|). The calculations do not reproduce very 
well the observed D^/D^ ratio and thus two different up-scale factors K 
would be needed to reproduce charged and neutral D cross sections simul- 
taneously. (See also a comment concerning this situation in Q). Since the 
spectral distributions for and {D~ and D ) are similar (see below), 
we use in the following an average K factor of 5.0 determined to reproduce 
the overall D yield. 

The hadro-production of open charm mesons on nuclear targets has re- 
cently been precisely measured in collisions induced by 800 GeV protons. 
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Figure 1: Cross section for forward (xf>0) D*^, D (left) and 1)+, D (right) 
production in proton and pion induced reactions, as a function of the center 
of mass energy. The data were taken from NA16 [|, NA27 @, E743 
E653 0, E789 0, E769 0, NA32 ||, WA92 ||], and E706 0. The 
sohd hues represent the pp Pythia calculations, up-scaled by a factor of 3.5 
(left panel) and 7 (right panel). 



Parameterizing the nuclear target dependence in the usual way, a{pA 
D + X) = (Ttv^", with (Tiv the open charm cross section for proton-nucleon 
collisions, the E789 experiment has obtained a = 1.02 it 0.03 it 0.02, for 
0.0 <xf< 0.08 (average xp of 0.031), comparing Be and Au targets |13]. 
With 250 GeV tt~^ and vr^ beams, the E769 experiment has measured 
a = 1.00 lb 0.05 lb 0.02 for D mesons with xp between 0.0 and 0.8, us- 
ing a target system composed of Be, Al, Cu and W foils |18]. Besides, E769 
has not found any significant dependence of a on or xp, at least for the 
a;F< 0.5 region, responsible for most of the cross section. These findings 
confirm what one expects for a hard process, in which the partons in the 
nuclear target independently contribute to the total cross section. 

In the following we will use hadron-hadron and hadron-nucleus data, 
along with the atomic mass number dependence, to evaluate to which extent 
charm production is described by leading order QCD calculations. To this 
end we have performed calculations with the code Pythia using MRS G 
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parton distributions |19| and {k^)=l GeV^/c^. We checked that using 
{k^)=2 GeV^/c^ changes the dilepton yield by less than 10%. Most of the 
calculations were performed for proton-proton collisions. We have, however, 
checked that isospin dependences are small for open charm production. At 
^/s = 18 GeV the differences predicted by Pythia in D-meson distributions 
from proton-proton and neutron- neutron collisions are below 10 %. Since 
open charm production is mostly due to gluon fusion this is not unexpected. 
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Figure 2: D+,D-,D°,D ,D+,D- production in p-p at 250 GeV: com- 
parison between experiment |^ (black dots) and the up-scaled Pythia 
calculation (histograms) . 



The energy dependence of charm production is rather well reproduced 
by these calculations (Fig. |T|). In Fig. 2 we present a comparison of Pythia 
predictions for pT and xp distributions with data from the E769 collabo- 
ration at a beam energy of 250 GeV |2C]. The absolute normalization of 
the calculations has been adjusted to reproduce the D meson data {D~^, 
D~ , D^, D , Df, D~), yielding an average K factor of 5.0, which is used 
for all subsequent calculations, independent of beam momentum. Notice 
that the value of K depends on the c-quark mass used; here we have used 
rric = 1.35 GeV/c^. It also depends on the choice of the scale used in 
Pythia. In our calculations we have used =s, the parton-parton center 
of mass energy squared. We have verified that, if the scale is given by 
rfij, (Pythia's default), the calculated cross sections increase by a factor 
of two, implying a corresponding decrease in the K factor. Here rfiT is the 
transverse mass of the outgoing particles produced in the hard scattering. 
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Note that the kinematical distributions do not depend appreciably on the 
choice of the scale. 

Apart from the absolute normalization, the Pythia calculations repro- 
duce reasonably well the measured and xp dependence, implying that 
the basic kinematics of the charm production process is correctly simulated. 

The agreement with the data points can be further improved by using 
a (^-function for the fragmentation of the charmed quarks in the charmed 
mesons. However, we find it difficult to believe that a J-function can be 
a good description of the hadronization process and, therefore, have used. 
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Figure 3: DD data from vr^-Cu collisions at y^=26 GeV [16|. Parts a) 
and b) represent the and xp distributions of the DD pair. Parts c) and 
d) show the DD invariant mass and the rapidity difference Ay = yo — Ujy 
distributions. The data are compared to the Pythia calculation (solid line). 
The Pythia distributions are normalized in the same way as the data, i.e. 
such that the integral over the displayed range is unity. 
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in the calculations reported in this paper, the Lund string fragmentation 
model (Pythia's default). 

To further check the kinematics of the charm production process as im- 
plemented in Pythia we present in Fig. 3 comparisons with the WA92 data 
for 7r~- Cu interactions, where kinematical correlations between charmed 
particles were measured at -^5=26 GeV \1G\. Fig. 3a shows good agreement 
between the transverse momentum spectrum calculated with Pythia and 
the data. As shown in Figs. 3b, 3c, and 3d, also the Feynman x distribution 
of the DD pair, the distribution in DD invariant mass Mj-,jy, as well as 



the rapidity difference between Z)+ (D^) and (D^) are reasonably well 
reproduced by the Pythia calculations, lending considerable confidence to 
the accuracy of the calculated lepton pair mass spectra. 
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Figure 4: Left: Invariant mass spectra of dielectrons from charm decay 
in pp collisions for different beam momenta (Pythia). Right: The same 
normalized to the spectrum at 450 GeV. 

To complete the systematics, and to lead into the discussion of charm 
contributions to dilepton spectra, we present in Fig. 4 the di-electron mass 
distributions from semi-leptonic D-decays, without any phase space selec- 
tion, for various beam energies. The calculated mass distributions are all 
rather broad, with a shallow maximum near M^+^- = 0.5 GeV/c^. No 
strong energy dependence is observed. To emphasize this point we show, 
in the right hand part of Fig. 4, the ratio of mass distributions relative to 
the one computed at 450 GeV. From this plot one observes only a slight 
narrowing of the mass distributions with decreasing beam momentum. Es- 
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pecially in the mass range < M^+^- < 1 GeV/c^, the variation of the 
ratio with mass is less than 20 %. 

The shape of the lepton pair mass distributions depends crucially on the 
implementation in Pythia of semi-leptonic D decays. In order to check 
if such decays are treated correctly we compare, in Fig. ^, the momentum 
distribution of electrons from D decays in Pythia with experimental distri- 
butions |21|. The theoretical distributions, where modes with K mesons and 
with K* mesons contribute with roughly equal weight, are in good agree- 
ment with the measured spectrum. Furthermore, we have verified that the 
values used in Pythia for the semileptonic decay branching ratios are in 
excellent agreement with the most recent values of Q. We conclude, there- 
fore, that uncertainties in the lepton energy spectra and branching ratios 
from D meson decays can be neglected in the following discussion. 
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Figure 5: Momentum distribution of electrons from and decays, 
measured in the rest frame of the decaying mesons. The experimental data 



(full dots) were taken from ref. [21|. The solid histogram was extracted from 
the Pythia calculation. 

The shape of the dielectron mass spectrum also depends on the shape of 
the DD invariant mass spectrum. In 1983, Fischer and Geist ||2^ proposed a 
scaling curve for the ratio between the width of the parent mass distribution 
and the threshold mass m*^'' = 2mD, plotted against the ratio between 
threshold mass and ^/s. Here, rriD is the mass of the D meson. This curve 



was motivated by ISR data and by early QCD calculations [23|. In Fig. 6 



we compare it with the present Pythia calculations. Although a trend 
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Figure 6: Relation between the collision energy and the half width at half 
maximum of the M^-^j distribution, as proposed by Fischer and Geist (solid 
line) and as calculated by Pythia for pp (stars connected by dashed line). 
Black dot: WA92 measurement [16|. 



similar to the Fischer-Geist curve emerges from the present calculations, 
the detailed shape looks quite different. We note, however, that the WA92 
measurement (black dot in Fig. 6) is in much better agreement with the 
present calculations than with the Fischer-Geist scaling. 



3 Dilepton yield from associated charm produc- 
tion in proton induced collisions 

In the previous section we have shown that calculations with the code 
Pythia describe reasonably well the most recent data onxp and px distribu- 
tions of D mesons as well as rapidity correlations and invariant mass spectra 
of the DD pairs [p^p^ in the center-of-mass energy regime of ^/s ~ 20 GeV. 
In addition, the characteristics of the semi-leptonic decays of D mesons are 
correctly reproduced. This gives us confidence that the level of lepton pair 
yield from charm production can be reliably calculated for pp and pA col- 
lisions on the basis of a Pythia simulation. We will now compare our 
estimate of the lepton pair yield from charm production to data on the in- 
clusive lepton pair yield. Both leptons from D meson and from Ac decays 
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are included in the calculations. 

First we discuss the results of Chilingarov et al. who have measured 
e^e~ pairs in pp collisions at the ISR {^/s ~ 60 GeV). For this energy 
Pythia (including the above discussed K factor) predicts a cc cross section 
of 100 /ib, roughly in agreement with the value discussed in [24|. The e'^e~ 
data of Chilingarov et al. and our calculation for e^e~ pairs from open 
charm are shown in Fig. 0. We have applied the cut > 1-6 + 0.5 M^+^- 
(with and M^+^- in GeV) on the electron pair transverse momentum in 
order to simulate the experimental cuts (see Fig. 5 in 




(GeV/c') 

Figure 7: Invariant mass distribution for e'^e~ pairs from pp collisions at 
■y/i=60 GeV, measured at mid-rapidity. The solid line represents the e'^e~ 
pairs from charm decays, as predicted by our Pythia calculation. 



According to the present calculations, and contrary to the conclusions 
of [24 1, the contribution from charm decays is approximately one order of 
magnitude below the observed dilepton yield while the shape of our com- 
puted e^e~ mass spectrum is rather similar to that displayed in Fig. 6 of 
ref. [24|. In order to check whether this result depends on the value of (A;^) 
used in the simulation, we repeated the procedure with (/c^)=2 GeV^/c^. 
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The corresponding yield was higher by only about 30%. Of similar size 
are the uncertainties due to the (small) differences between calculated and 
measured xp distributions for D mesons (see Fig. 2). We assume that the 
difference between our results and the conclusions of Chilingarov et al is con- 
nected with the fact that, in ref. [24|, the yield for e'^e~ pairs was calculated 
assuming isotropic decay in the pair rest frame, while Pythia incorporates 
a strong anisotropy. In addition, in ref. [^^, it was assumed that the lepton 
pairs are produced with uniform rapidity distribution. The fact that Pythia 
reproduces well the observed distribution of rapidity gaps between the DD 
mesons (see Fig. 3), which is at the root of the anisotropy, lends strong sup- 
port to the validity of the present calculations. We therefore conclude that, 
at ISR energies, the charm decay contribution to the dilepton continuum 
below the (/)-mass is small. To explain the observed continuum one should 
probably revisit the Dalitz decay contributions to the mass spectra. 




Figure 8: Left: Invariant mass spectra of dielectrons from charm decays in 
p-p collisions at beam energy 450 GeV (Pythia), before (upper histogram) 
and after (lower histogram) applying the CERES filter. Right: the ratio of 
the two histograms. 

Recently the CERES/NA45 experiment has measured e^e~ pairs in p-Be 
and p-Au collisions at 450 GeV/c incident proton momentum [^. In order 
to compare our Pythia calculation to the CERES data, we included in our 
simulation their momentum resolution and kinematic cuts (pT > 50 MeV/c 
and 2.1 < rj < 2.65 for electrons and positrons). The effect of the CERES 
acceptance on the e~^e~ mass spectra from charm decays is severe, as can 
be seen from Fig. |8|, where we have compared e'^e~ mass distributions from 
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charm before and after applying the acceptance cut. The integral of the mass 
distribution is reduced by approximately a factor of 60 and, in addition, 
the shape is distorted, leading to a steeper mass dependence for large pair 
masses. The dependence of the cross section on the nuclear mass number was 
taken into account as outlined in the next section. Finally the cross section 
was converted into a multiplicity distribution and scaled by < dNch/dy 
taken from the CERES experiment which gives the same normalization as 
used by the CERES collaboration. 

Fig. 1^ shows the CERES data as published in ref. |^], but including 
our estimate of the contribution from charm production. Clearly, also in p- 
nucleus collisions charm production does not give an important contribution 
to the dilepton continuum below the light vector mesons. 
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Figure 9: Invariant e~^e~ mass distribution from 450 GeV/c p-Be and p-Au 
collisions measured by the CERES experiment. The background from asso- 
ciated charm production has been added to the hadron decay contribution 
according to our estimate. 

HELIOS/NA34-1 |25| has published a measurement of e~^e~ and fJ-^fJ^^ 
pairs in p-Be collisions at 450 GeV/c momentum for the same mass region. 
From their analysis one concludes that the charm decay contribution to their 
data is negligible. 
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While lepton pairs from associated charm production can be neglected 
for low invariant masses, for masses above the i;^i-meson the situation is 
quite different. Two experiments - NA38 and HELIOS/NA34-3 - have 
measured muon pairs from p-nucleus collisions in this mass range. Both 
experiments have shown that the dimuon cross section in p-W collisions at 
200 GeV/c can be accounted for by charm production and pairs produced 
via the Drell-Yan mechanism. The HELIOS-3 estimate of the charm 
cross section was based on an approach similar to the one presented in 
this paper. The experiment NA38 used Pythia to determine the shape 
of the dilepton background from charm production but fitted the charm 
cross section to their muon pair data [^. The result is compatible with the 
cross sections shown in Fig. |l[ Both experiments consistently show that, for 
masses above the (j) meson, cc-production is an important source of lepton 
pairs until the continuum is overwhelmed by pair production via the Drell- 
Yan mechanism in the mass region beyond the J/V'. 



4 Extrapolating to nucleus-nucleus collisions 

Since the charm cross section scales with the number of target nucleons going 
from p-p to p-nucleus collisions, we scale the p-p calculations from Pythia 
with the total number of nucleon-nucleon collisions to obtain the nucleus- 
nucleus cross section. As discussed above, isospin effects can be neglected 
for the calculation of open charm production within an accuracy of better 
than 10 %. The total number of nucleon-nucleon collisions is calculated from 
the nuclear collision geometry based on the concept of an overlap function 



TAsib) ioT A + B collisions at impact parameter 6 |27|. This leads to: 



dy dy 



/ d%TAB{b). (1) 
Jbi 



Here, the integration limits indicate the impact parameter range covered in 
the experiment. For minimum bias collisions, where hi = and b2 = oo this 
leads to 

^■'"=AB^" (2) 

dy dy 

For the rapidity density distributions this implies 

dy dy '^^"'^^ vr(62-62) • W 
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Here, crf^^^ is the inelastic nucleon-nucleon cross section at the appropriate 
energy 0. Equation |3|, when apphed to central collisions (6 ~ 0) and using a 
sharp sphere approximation of the nuclear density profile for the evaluation 
of the overlap function, results in the familiar dependence 

^ (6^0)«^ A^/^. (4) 
dy dy 

We note that this approach neglects possible effects due to the Fermi motion 
of the nucleons in the colliding nuclei. However, the threshold for open charm 
production is -y/i = 5.5 GeV compared to the -y/i = 17.3 GeV available in 
Pb-Au collisions so that Fermi motion effects are not expected to play an 
important role. We use Woods-Saxon distributions for the nuclear densities 
to evaluate the thickness function TAB(b) (or Ta for the case of pA 
collisions) and adjust the impact parameter range to obtain the fraction of 
inelastic cross section sampled in each experiment. Applying this procedure 
to the pp results from Pythia we obtain the D meson yield expected for 
nucleus-nucleus collisions. In Fig. |l^ we show the prediction for the rapidity 
distribution of D mesons in central (5 % of geometrical cross section) Pb+Pb 
collisions at 160 GeV/c/nucleon. The distributions exhibit rather drastic 
differences between D mesons and their anti-particles. 

These differences originate, within the Pythia framework, from the in- 
teraction of c and c quarks with valence quarks of the colliding nucleons 
which takes place prior to hadronization. While the c forms a string with a 
single valence quark, the c interacts with valence diquarks, and is dragged 
towards beam or target rapidity. As a consequence, the and meson 
rapidity distributions are broadened and, at energies around -y/s =20 GeV, 
even exhibit minima at midrapidity. 

Similar asymmetries, although different in detail because of the differ- 
ent quark content of the projectile, have recently been observed for tt~ - 
nucleus interactions at 500 GeV/c [29|. There it is also shown how well 
these asymmetries are reproduced by Pythia^ It would be interesting to 
verify experimentally whether such asymmetries also exist in proton induced 
collisions, in particular at y/s ~ 20 GeV, since according to Pythia these 
asymmetries tend to disappear at higher beam energies. 



^ In actual calculations the inelastic nucleon-nucleon cross section drops out as we 
compute dir^ly^- = If- a-, 

^ In ref. [M it is shown that the measured asymmetries can be well reproduced by 
Pythia if one sets the mass of the charmed quark to 1.7 GeV/c^ and otherwise uses the 
same parameters as are used in the present investigation. With our parameter set the 
agreement between data and calculations is only slightly deteriorated. 
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Figure 10: Calculated rapidity distributions oi D^, , D^, from central 
(5%) Pb+Pb collisions at 160 GeV/c/nucleon. The ordinate represents 
absolute multiplicities per unit rapidity. 



The shapes of the D meson rapidity distributions for Pb+Pb collisions 
might also be distorted in the nuclear environment, as observed for dis- 
tributions measured in p- nucleus collisions, commonly referred to as Cronin 
effect [30j. Clearly, a direct measurement of D meson production in Pb+Pb 
collisions would be of great interest. The possible modifications discussed 
above should not change the total rate significantly and we expect that (in 
the absence of charm enhancement) the rapidity density of D mesons from 
Pb+Pb collisions at 160 GeV/nucleon is dA^/dy=0.10-0.15 at midrapidity. 

As in the previous section we want to focus on the charm decay contri- 
bution to the inclusive lepton pair production. In Fig. 11 the results of our 
calculation are compared to the S-Au e~^e~ data published by the CERES 
group p . Clearly, charm production in the "standard" scenario as presented 
here contributes less than one percent of the observed dilepton yield. The 
slight difference in shape of the charm contribution compared to Fig. |9| is 
solely due to the higher px-cut used by the CERES experiment for S-Au 
collisions. 

It is not unreasonable |^l] to contemplate the possibility of a strong ad- 
ditional source of charm in nucleus-nucleus collisions. However, an enhance- 
ment by a factor of 150 would be needed to describe the present CERES 
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Figure 11: Invariant mass spectrum of e^e^ pairs measured by CERES 
in 200 GeV/nucleon S-Au collisions. The solid line gives the contribution 
expected from charm production in the standard scenario. The dashed line 
results from adding to the hadron decay background a charm contribution 
enhanced by a factor of 150. 




Figure 12: Dimuon spectra from p-U, S-U, and Pb-Pb collisions 
(NA38/NA50 Collaboration). The S-U spectrum excludes charm enhance- 
ment by more than a factor of 3. The figure was taken from |32| and is 
based on results from [Esl- 
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data. The result of such an enhancement is shown as the dashed hue in 



Fig. 11. Though it roughly explains the overall dilepton yield, it predicts 
a flatter mass distribution than is observed resulting in an overestimate of 
the yield at high masses. The results for dimuon production obtained by 
the NA38/NA50 Collaboration are shown in Fig. |l^ where also the predic- 
tions for the charm and Drell-Yan continuum are indicated. Based on these 
results it is clear that charm enhancement factors larger than 3 would lead 
to inconsistencies with the S-U dimuon data. 



5 Conclusions 

In summary, the level of charm production in p-p and p-A collisions in the 
energy range ^/s ~ 20 — 40 GeV is now reasonably well measured. The 
trends in the data are well reproduced by leading order perturbative QCD 
calculations performed using the code Pythia with {k^)=l GeV^/c^. The 
calculations were extrapolated to nucleus-nucleus collisions by scaling with 
the geometrically allowed number of nucleon-nucleon collisions. We have 
used these calculations to determine the dilepton pair contribution expected 
from associated cc production. In the mass region below the vector mesons 
p, uj and (f) this contribution is negligible compared to the yield arising from 
decays of light mesons for all collision systems investigated in this paper. 
For S-Au collisions the prediction is two orders of magnitude below the 
experimentally observed e"^e~ continuum. Since any possible enhancement 
of charm production in nuclear collisions is limited by p.'^ p~ data in the 
mass region above the cj) to less than a factor of 3, any explanation of the 
observed anomalous low mass pair production based on charm production 
can be ruled out. 

We thank George Ginther, Karl Harrison and Thomas Ullrich for use- 
ful discussions while collecting the charm cross section measurements from 
different experiments. We are grateful to H. G. Fischer for helpful and 
provocative remarks about charm production and to W. Geist, M. Mangano 
and A. Morsch for helpful discussions. 
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